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Ferromagnetic, in contrast to antiferromagnetic, Heisenberg chains can undergo a Spin-Peierls 
dimerization only at finite temperatures. They show reentrant behavior as a function of temperature, 
which might play a role for systems with small effective elastic constants as, for example, monatomic 
chains on surfaces. We investigate the physical properties of the dimerized ferromagnetic Heisenberg 
chain using a modified spin-wave theory. We calculate the exponentially decaying spin and dimer 
correlation functions, analyze the temperature dependence of the corresponding coherence lengths, 
the susceptibility, as well as the static and dynamic spin structure factor. By comparing with 
numerical data obtained by the density-matrix renormalization group applied to transfer matrices, 
we find that the modified spin wave theory yields excellent results for all these quantities for a wide 
range of dimerizations and temperatures. 
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I. INTRODUCTION 

One-dimensional (ID) spin systems play an important 
role in quantum magnetism. A basic model is the an- 
tiferromagnetic (AF) spin- 1/2 Heisenberg chain which is 
exactly solvable by Bcthe ansata 1 - and has gaplcss excita- 
tions (spinons). Many realizations of quasi ID AF spin- 
1/2 Heisenberg chains are known today with Sr2Cu0 3 
being one of the best studied examples^ An analysis in 
terms of a spin-only model, however, might not always 
be applicable because of a coupling to lattice degrees of 
freedom. In complete analogy to the well-known Peierls 
effect 3 - for ID conductors^— a spin-Peierls transition can 
occur leading to a dimerization of the spin exchange and 
a gapped excitation spectrum*^ A well known example 
is the spin-Peierls transition in CuGcOsi^ 1 - 

In contrast, ferromagnetic (FM) spin chains are less 
frequent. An example which has been analyzed in some 
detail is the organic system (CH^NCuCb}! 1 ^ The mag- 
netization and susceptibility data of this system are well 
described by the spin-1/2 FM Heisenberg chain. The 
thermodynamic properties of this fundamental model can 
be calculated by Bethe ansat a 13 : 14 as well as by a mod- 
ified spin- wave theory (MSWT)* 1 ^ Further examples for 
realizations of a ID Heisenberg model with FM exchange 
are CuCl 2 (DMSO) and CuCl 2 (TMSO)^ where DMSO 
(TMSO) stands for dimethylsulfoxide (tetramethylsul- 
foxide), (C 6 HilNH 3 )CuCl 3 and (C 6 H 1 lNH 3 )CuBr 3? i2as 
well as 2-benzimidazolyl nitronyl nitroxide (2-BIMNN).— 

Another fascinating system are ID arrays of Co chains 
which sclf-asscmblc on Pt substrates. The magnetic ex- 
change in these nanostructured monatomic chains has 
been shown to be predominantly FM . 19 ' 20 The spe- 
cific spin texture in such systems, however, might be 
complicated^ due to the reduced symmetry allowing for 
exchange interactions different from a pure Heisenberg 
exchange. Since the atoms can be easily moved along the 



surface, the effective elastic constants of these monatomic 
chains are expected to be small. This has lead us, inter 
alia, to the question if a coupling to lattice degrees of 
freedom can induce a dimerization of FM exchange in- 
teractions similar to the AF case. 

In a recent work,— we have shown that such a spin- 
Peierls effect is indeed possible for the FM Heisen- 
berg chain — at least at the level of the adiabatic 
approximation — but has to be activated by thermal fluc- 
tuations. Furthermore, we have argued that this mech- 
anism seems to explain the observed dimerization of the 
FM exchange in the C-phase of the perovskite YV0 3 ) 23 i 24 
with the dimerization being caused by a coupling to or- 
bital rather than lattice degrees of freedom. 

The aim of this paper is to study the physical prop- 
erties of a dimerized FM spin chain independent of the 
mechanism which causes the dimerization. As an an- 
alytical method to calculate the thermodynamic prop- 
erties of the dimerized ferromagnetic Heisenberg chain 
we will use a MSWT and compare with numerical data 
obtained by the transfer matrix renormalization group 
(TMRG). As a central result of our paper we show that 
the MSWT yields excellent results for a wide range of 
dimerization parameters and temperatures, i.e., when we 
compare with TMRG data. In particular, we find that 
at finite temperatures T both the spin and the dimer 
correlation function decay exponentially with the same 
correlation length (we set h = fee = 1 ) 



£ s = c(T,6)JS 2 /T, 



(1) 



where J < is the FM exchange constant in the Heisen- 
berg Hamiltonian and S is the spin quantum number. 
The lattice spacing has been set to unity. A detailed dis- 
cussion of the function c(T, S) which captures the effects 
of the dimerization S as well as the higher order temper- 
ature dependence of the correlation length is given for 
5 = 1/2 and S = 1. Moreover, we also present explicit 



analytical expressions for the doping and temperature de- 
pendence of the asymptotic behavior of correlation func- 
tions and the spin susceptibility. 

Then we address the question how the spin-Peierls 
symmetry breaking manifests itself in the spin-structure 
factor given that the correlation functions decay expo- 
nentially. To detect this symmetry breaking we show 
that it is useful to study an off-diagonal spin-structure 
factor. Three approaches are used to calculate the static 
structure factors: (i) via the dynamic spin structure fac- 
tors, and via the equal-time spin correlation functions as 
obtained by (ii) MSWT and (hi) TMRG. 

The paper is organized as follows: In Sec. [Ill we intro- 
duce and motivate the model. Furthermore, we discuss 
the set-up of a MSWT which will be used throughout 
the paper to obtain analytical results for various thermo- 
dynamic quantities. In Sec. Hill spin correlation functions 
for the uniform and the dimerized ID chain are discussed. 
Here we derive analytic expressions for the correlation 
functions both in the short and the long-range limit, cal- 
culate the magnetic susceptibility and compare with data 
obtained from numerical TMRG calculations. In Sec. IIVI 
and Sec. [V] we discuss the dynamic and the static spin 
structure factor, respectively. In Sec. IVII we present a 
short discussion and a summary of our results. 



II. SPIN-PEIERLS PHASE IN ID 
FERROMAGNETS 

The Spin-Peierls Hamiltonian we want to investigate 
reads 



iV 



H 



(2) 



J^{l + (-l)^}S r S i+1 , 
i=i 
where J > is the nearest neighbor exchange interaction 
among the spin S operators Sj and Sj+i on the sites j 
and j + 1, respectively, and S <= [0, 1] is the dimcrization 
parameter. Consistent with dimerization we consider an 
even number TV of spins and assume periodic boundary 
conditions. In Ref. [22| the Hamiltonian, Eq. ([2]), was in- 
vestigated by TMRG. This method is particularly suited 
to study the properties of ID quantum systems directly 
in the thermodynamic limit. Details about this approach 
are given in Refs. [25H281 In particular the gain in the free 
energy 6f ~ —T 3 ' 2 5 2 stemming from a dimerization was 
compared in an adiabatic approximation to the cost in 
elastic energy ~ Kd 2 caused by the lattice distortion^ A 
dimerized phase was found to be stable in a regime of fi- 
nite temperatures for a sufficiently small clastic constant 
K < K c , with K c /J ~ 0.118. 

In order to apply the MSWT, we introduce boson op- 
erators {aj„, a m } on site m via a Dyson-Maleev transfor- 
mation 
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FIG. 1: (Color online) Magnon dispersions obtained from 
MSWT for various dimerizations. The lower branches cor- 
respond to LjjT while the upper branches are given by uj^ . 



Using the Dyson-Maleev transformation for the Hamil- 
tonian, Eq. ([2]), and expanding to bilinear order in the 
boson operators, H is easily diagonalized in terms of new 
boson operators {a kl ct k ; f3 k , f3 k } obtained by subse- 
quent Fourier and Bogoliubov transformation. The in- 
troduction of two sets of bosonic operators is required 
here since the unit cell of a dimerized FM chain is dou- 
bled. We find 
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Here Eg = —JNS 2 is the ground state energy at zero 
temperature (T = 0), which is independent of the dimer- 
ization (5, k is the ID momentum, and uj k stands for the 
dispersion of the two magnon branches given by 



2JS(l±Vcos 2 k + 5 2 ', 



'k 



(4) 



From this expression it becomes clear that a finite dimer- 
ization, 5 > 0, induces a splitting of the spin- wave dis- 
persion into an accoustic and an optical branch at the 
boundaries k = ±ir/2 of the Brillouin zone. As a conse- 
quence, the two branches flatten with increasing dimer- 
ization. This is shown in Fig. [1] where the dispersions for 
various values of <5 are shown. 

To fulfill the Mcrmin- Wagner theorem of vanishing 
magnetization at finite temperature T, usual spin-wave 
theory has to be modified by a Lagrange multiplier serv- 
ing as a chemical potential us (T) that sets the finite tem- 
perature magnetization to zeroi 15 i 30 This results in the 
constraint 
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Here 
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FIG. 2: (Color online) Occupation numbers for the (a) lower and (b) upper branch for S — 1 and T — 0.5J and various 
dimerizations as indicated. 



is the Bose distribution function, /3 = 1/T, and 



(7) 



is the reduced magnon dispersion. The chemical poten- 
tial fis(T) will be negative and vanishes as T — > 0. 

Results obtained from this procedure are in excellent 
agreement with the Bethe ansatz results for the 5 = 1/2 
uniform ferromagnetic (5 = 0) if ts/S < 1, where we 
have defined a reduced temperature 



frf 
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JS(l-<5 2 )■ 



(8) 



The occupation numbers for the lower and upper branch 
n B(C/7) an d n B(Cfc") are shown in Fig. [2] for T/J = 0.5 
and 5 = 1. As expected, the occupation for the lower 
branch shows a maximum at k = whereas the occupa- 
tion for the upper branch has a maximum at k — tt/2. 
With increasing dimcrization these peaks become less 
pronounced due to the flattening of the dispersions, see 

Fig. m 



III. CORRELATION FUNCTIONS 

To obtain the spin correlation function (SCF) for an 
arbitrary distance \i — j\ it is crucial to include also the 
quartic terms in the bosonic operators stemming from 
the Dyson-Malecv transformation of the scalar product 



Si • S,. This leads to 



(8i8j) 



5(5 + 1) iori=j 

\{[a\aj) + (atoi)J iori^j 



(9) 



For T > the correlation function approaches zero in the 
limit \i — j | — ^ oo as (a]aj) vanishes in this case. 
For the dimerized chain we obtain from Eq. ([9]) 
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FIG. 3: (Color online) Nearest neighbor spin correlation func- 
tions on (a) the strong bonds and (b) the weak bonds as a 
function of dimerization for various temperatures and 5=1. 
The solid lines show the MSWT results. For comparison, 
TMRG data are shown as symbols. The inset shows data for 
low temperatures. 
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FIG. 4: (Color online) Spin correlation function (solid lines) 
and alternation correlation function A(r) (dashed lines) for 
S = 1 and S — 0.1 and various temperatures as a function 
of distance r obtained by MSWT. For a given temperature, 
the exponential decay of both quantities is determined by 
the same correlation length. 
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FIG. 5: (Color online) Spin correlation function (solid lines) 
and alternation correlation function (dashed lines) obtained 
by MSWT for S = 1 and T/J = 0.1 and various dimeriza- 
tions as indicated. 



if r 7^ is even, and 

(Sj • S j+r ) = (^ E [»B(C fc -) - MO] fk(j,r,5) 
if r is odd. Here we have denned 



odd r as 



(11) 



, . „ cos(fcr) cosfc + (— lV(5sm(£;r)sinfc 

fk(3, r, S) = === . (12) 

Vcos 2 k + S 2 sin k 

Note that the summation in Eqs. (|1Q|) and (|11|) runs over 
the reduced Brillouin zone which is the reason for the 
minus sign in Eq. (|11[) . From these expressions we ob- 
tain the nearest neighbor SCF shown in Fig. [3] One 
finds that the correlations on the strong bonds are en- 
hanced while those on the weak bonds are reduced when 
the dimerization is increased. The dimerization strongly 
affects the SCF on the weak bonds whereas the strong 
bonds arc far less affected. We find excellent agreement 
between the results obtained by MSWT and TMRG data 
up to T/J = 1 and even the value at full dimerization 
(5 = 1) is captured correctly. However, MSWT suggests 
a quadratic behavior in e = 1 — S on the weak bonds for 
e <C 1 whereas a perturbativc treatment starting from de- 
coupled dimers shows a linear behavior as confirmed by 
the TMRG data. Finally the steep decrease of the nearest 
neighbor SCF on the weak bonds for low temperatures 
and strong dimerization indicates a non-analyticity for 
T = at 5 = 1. Thus the ground state of the dimerized 
chain is always a uniform ferromagnetic state except for 
S = 1 where it consists of decoupled dimers. 

Next we shall explore the long-distance behavior of 
the spin-correlation function (Sj ■ Sj +r ) and its varia- 
tion with dimerization S and temperature. In addition 
we define an alternation correlation function (ACF) for 



A(r) 



= I (Sj ■ Sj+r, 



(Sj ■ Sj 



(13) 



This expression describes the alternation of the SCF be- 
tween the weak and the strong bonds and is thus a mea- 
sure for the dimerization of the system. We exclude even 
r in Eq. (|13|). as in this case A(r) = by symmetry. In 
Fig. H the SCF and the ACF obtained from MSWT for 
5 = 1 and S = 0.1 are shown for different temperatures 
and in Fig. | for S = 1 and T/J = 0.1 and different 
dimerizations. The asymptotic behavior of the correla- 
tion functions is governed in both cases by an exponential 
decay ~ e~' r !/^ with the same correlation length £$■ As 
expected, increasing dimerization 5 reduces the correla- 
tion length. For large dimerization the alternation of the 
exchange integrals leads to a staircase-like behavior of 
the SCF (Sj ■ S j+r ) (see Fig. E for 6 = 0.6). 

The function c(T, 6) (see Eq. (JT])) describes the behav- 
ior of the correlation length as a function of dimerization 
and also gives corrections to the leading order 1/T tem- 
perature dependence. For a uniform ID ferromagnet to 
lowest order in T we have c(T, 5 = 0) = 1j22 In order to 
go beyond this limit we have to investigate in detail the 
behavior of the SCF at large distances. For sufficiently 
low temperatures we may neglect contributions stemming 
from ^b(C^) in Eqs. ©, (|TU)) . and (fTTj) . Expanding w^T 
and fk (J, r, S) to lowest non- vanishing order and using a 
saddle point integration in Eqs. (fT0]l and (fTTj) we obtain 
in a large-r expansion (r » 1) 
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TABLE I: Coefficients of c(T, S) given in Eq. {THJ as extracted 
from Eq. ifl?) for S = 1/2 and S = 1. The coefficient oi is 
given by (1 — S 2 ). 



with the correlation length 
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Here v$ = —[is/T is the negative reduced chemical po- 
tential reading 



ts ( i — 3 2 a3 3/2 



, (16) 



with a = C(n )/(^V^) where ((z) is Riemann's zeta func- 
tion. Eq. ([It]) explains why the SCF and the ACF both 
decay with the same correlation length. We note that for 
<5 = this expression agrees with the one obtained for 
the uniform case^i 

Inserting Eq. (|16|) into Eq. (fl5|) and expanding for tg <C 
1 we obtain the analytical expression 



c(T,d)^(l-5 2 ) 
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In Fig. [SI c(T, 5) for S = 1 and various dimerizations 
is shown in comparison to numerical data from TMRG. 
We find excellent agreement between both methods. In 
Fig. c{T, 5) is shown for S = 1/2 and S = 1 as a 
function of temperature for 5 = 0.1. Note that c(T,S) is 
much steeper for S = 1/2, sec Eq. ([17]). In Tab. [J the 
coefficients ai, 02, 03, and 04 of the scries 



c(T, <5) w ai + a 2 \/ — + a 3 



"4 



5/2 



(18) 



are shown as extracted from Eq. (flTj) . We remark that 
the higher order corrections 03 and 04 are much larger 
for S — 1/2 than for S 1 = 1 and increase with increasing 
dimerization. 

Apart from the large-r expansion, an analytical ex- 
pression for small distances can also be obtained from 
Eqs. (fT0|) and (fTTj) . Again neglecting contributions stem- 
ming from n,B(Ck ) ano - expanding the dispersion as well 
as the function fh{j, r, 6), Eq. (fi"2|). to second order in k 
we find the short-distance expansion 
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/J (19) 

In Fig.[5]the SCF is presented for S = 1, T/J = 0.1 and 
<5 = 0.3 as a function of distance. Fig. [8ja) compares a 
numerical self-consistent solution of Eqs. (flO|) and (fTTj) 
with the expansion for small distances r, Eq. (|19|) . and 
the long distance asymptotics, Eq. (fT4")l . Furthermore, 
numerical TMRG data are shown. For r S> 1 the ex- 
ponential decay of the SCF, as derived from the MSWT 
is clearly visible. When applicable, we find good agree- 
ment between the analytical expressions for short and 
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FIG. 8: (Color online) (a) Spin correlation function 
{Sj ■ Sj+r) on a logarithmic scale for S — 1, T/J = 0.1 
and 8 = 0.3 obtained by TMRG compared with the ana- 
lytical results from the asymptotic expansion, Eq. (|14|) . and 
the short-range expression, Eq. (|19|) . up to r — 4. The solid 
line connects data obtained from a numerical solution of the 
MSWT, Eqs. (Unjl and (TJT]), for integer r. (b) Same as above 
on a linear scale. 
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FIG. 9: (Color online) Difference between the spin correla- 
tions on weak and strong bonds, see Eq. (|13[) . for S — 1, 
T/J = 0.1, and two representative dimerization parameters. 
The filled symbols are data obtained from TMRG whereas 
the open symbols show results obtained by MSWT. 
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FIG. 10: (Color online) Inverse susceptibility 1/x as a func- 
tion of T/J for S — 1 and different values of 5 as indicated. 
The lines show the results from MSWT, the symbols for 
8 < 1 TMRG data, and the symbols for 8=1 the exact 
solution for decoupled dimers. An increasing dimerization 
suppresses \, and the decoupled dimers follow a Curie law 
X ~ 1/T. 



long distances and the results obtained by TMRG. The 
same quantities for the same parameters are shown in 
Fig-E^b) on a linear scale. In Fig.[9]wc compare the ACF 
for spin 5=1 and two large values of the dimerization 
parameter, S = 0.6 and 0.9, as obtained from MSWT, 
with TMRG data. Although a high dimerization leads 
to a high effective temperature t$ the agreement of the 
MSWT with the numerical data is still remarkably good. 



Next, we calculate the zero field susceptibility per site 

WrY.&rSx). (20) 
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Applying the Dyson-Maleev transformation we find 
S 
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In order to obtain an analytical result for the susceptibil- 
ity, we may interpolate the SCF in Eq. (|20|) between the 
short-range limit given in Eq. (|19p and the asymptotic 
expression of Eq. (|14[) . From this we obtain 
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(22) 

where we have used that to lowest order the correlation 
length is given by £5 = S/tg. In the limit of vanishing 
dimerization our result reduces to the result previously 
found by Takahashi for the uniform chain i 30 ' 32 

In the opposite limit of full dimerization the suscepti- 
bility can easily be calculated exactly. In this case we 
have decoupled ferromagnetic dimers and the sucepti- 
bility follows a Curie law at low temperatures. Within 
MSWT we find from Eq. fl2l]) 



X(S = 1) 



S + 2{nl(v 1 ) + n 2 s (AJS + v 1 )} 
3T 



For T I J -^ 1 this expression reduces to 

S(2S +1) 



3T 



(23) 



(24) 



which is the Curie law (per spin) for decoupled ferromag- 
netic dimers. In Fig. [10] the inverse susceptibility 1/x as 
obtained from MSWT is compared with TMRG data. 
The agreement is good for ts/S < 1 and also the case of 
decoupled dimers is accurately captured by MSWT. 



IV. DYNAMICAL SPIN STRUCTURE FACTOR 

For the dimerized AF Heiscnbcrg chain the dynam- 
ical spin structure factor (dynamic SSF) has been in- 
tensely studied by various methods.— ~— Here we will 
discuss the impact of a finite dimerization on the dy- 
namic SSF of the FM Heiscnbcrg chain. The dynamical 
properties of the uniform ID ferromagnet within MSWT 
have been adressed by TakahashL 37 ' 38 In this case for a 
given value of q the dynamic SSF exhibits a magnon peak 
from which the magnon dispersion can be obtained. In 
this approximation magnon-magnon interactions are not 
taken into account and the peaks are broadened due to 
thermal fluctuations only. The magnon peaks, which ex- 
hibit a Lorenzian lineshape at finite temperatures, thus 
reduce to ^-functions at T = 0. Moreover, there is a 
two-magnon continuum beyond which the dynamic SSF 
vanishes. New aspects to this picture have been added 
recently 3 ^ when it was shown that an edge singularity 
occurs at the boundary of the two-magnon continuum 
caused by a diverging density of states. 

Here we want to generalize this analysis to the dimer- 
ized FM chain. Within the standard approach of neutron 



scattering the differential cross-section is directly related 
to the dynamic SSF42 The latter can be obtained by a 
Fourier transform of the two-point correlation function. 
It was shown in the context of X-ray scattering that fur- 
ther insight into the microscopic properties of physical 
systems may be achieved if the incident photon fulfills 
the Bragg condition such that the resulting state is given 
by a coherent superposition of the forward diffracted and 
the Bragg reflected wave4i~— The structure factor can 
be written as a matrix in terms of the reciprocal lat- 
tice vectors in analogy to the dielectric matrix whose in- 
verse defines the dynamic charge structure factor of a 
dielectric 4i~— 

In the same spirit we formulate the dynamic SSF of 
the dimerized FM chain as a 2 x 2-matrix 



S(q,ui) 



S(q,q,uj) S(q,q + tt,uj) 

S(q + n, q,w) S(q + n, q + ir, u) 



In order to obtain this matrix we calculate the Green's 
function 

G(j,r,r)^-(r[^(0)-5 J+r (r)]), (25) 

which within MSWT reads 

G(j,r,T)=-(T{a j (0)a] +r (T)}] 



T{ot(0)o J+r (r)} 



(26) 



From this the matrix elements of the dynamic SSF can 
be obtained straightforwardly via 



S(q,q',u>) = 2n B (-w)ImG rot (<3', q',u), 



(27) 



where G Iot (q,q' ,tu) is the retarded finite temperature 
Green's function obtained from G(j,r,r). For the di- 
agonal elements we find S(q,q 7 uj) = S(q + ir,q + n,u) 
with 

S{q,q,w)=^n B (-w) £ £ 5(u> - ef(k)) 

o\<t'S{±} A: 



x{l + aa'^(k,q + k)}Af{^,C +k ,p], 



(28) 



where we have defined 

cos k cos k' — S 2 sin k sin k 1 



$(k,k') = 



1 - (1 - S 2 ) sin 2 fcyT - (1 - S 2 ) sin 2 k' 

(29) 



and 



W, C +fc , 0\ ^MC k )n + MC +k )} 



( e -*r'«-i) 



(30) 



with the two particle continuum t a ° (fe) = Cfc — C+fc- ^ s 
in the case of the uniform ID ferromagnet the dynamic 
SSF for the dimerized chain fulfills detailed balance. 
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FIG. 11: (Color online) Diagonal elements of the dynamic 
spin structure factor S(q,q,Lu) for a dimerized ferromagnetic 
chain and S — 1, T/J = 0.1: (a) S(g, q, cu) spectra for different 
momenta q and 5 = 0.6. The red dashed lines show the edge 
boundaries and the symbols show the peak positions projected 
onto the (g,a>)-plane. (b) broadened spectra Ia{q, Q, w) for 
q = 4tt/5, 8 = 0.4 and for different A = 0,0.001,0.01, and 
0.1J. Vertical lines in A = panel indicate positions of edge 
singularities. 



In this approximation the two magnon continuum de- 
termines the regions where S(q, q, uS) is nonzero. This 
can be seen in Fig. [TTT a) where the diagonal element of 
the dynamic SSF for the dimerized chain is shown for 
8 = 0.6. The symbols reflect the peak positions pro- 
jected onto the (q, oj)-plane. They follow the reduced 
magnon dispersions £r. Since we are dealing with two 
separate magnon branches, an edge singularity occurs at 
each boundary of the two magnon continua due to the di- 




FIG. 12: (Color online) Off-diagonal element of the dynamical 
spin structure factor matrix S(q, q + 7r, w) for q = n/5, S — 
1, T/J = 0.1, and 8 as indicated. Note that the data for 
8 — 0.6 are shifted by 20 units along the ordinate for clarity 
of presentation. The vertical lines indicate positions of edge 
singularities (see text). 



vergence of the respective density of states. The bound- 
aries beyond which the dynamic SSF for the respective 
branch vanishes are given by the red dashed lines. For 
uj > wc find three boundaries out of the four possible 
combinations because for e~ + the argument of the delta 



function in Eq. (|28|) never vanishes. 

Experimentally, the energy resolution in a measure- 
ment is always finite. An important question to ask is 
therefore if the edge singularities can be resolved at all. 
We address this question by studying the dynamical SSF 
with a finite energy resolution A 

oj+A 

I A (q,q,u) = — / duj S(q,q,u ). 

ui— A 

Results for this quantity are shown in Fig. Illf b) for 
8 = 0.4 and various resolutions. Being very close to a 
magnon peak, the lower edge singularity has a substantial 
spectral weight. Even with a relatively moderate energy 
resolution this edge singularity is therefore detectable. 
On the contrary, the two edge singularities at higher en- 
ergies carry almost no spectral weight in the considered 
example making them experimentally irrelevant. 

The off-diagonal elements of the dynamic SSF, S(q, q+ 
tt,u>) = S(q + TT,q,co), are given by 

S^q + n^^nvi-u) ]T £ 5(u - ef(k)) 

<T,cr'e{±} k 

xaa'v&(fc,q + fc)AA[C fc CT ,C 9 CT + fc ,/3], 

(31) 
where we have defined 



#(fc,fc') = 



28 cos k sin k' 



yjl-{l- 8 2 ) sin 2 kJl - (1 - 8 2 ) sin 2 k' 



(32) 




From this we have 



FIG. 13: (Color online) Diagonal static spin structure factor 
for S = 1, T/J = 0.1. Symbols denote data from TMRG, 
lines results from MSWT. Main: Data for dimerizations 5 — 
(diamonds), 5 = 0.3 (triangles), S = 0.6 (squares), and S — 
0.9 (circles). Left inset: High-momentum peak for 5 — 0.9 
at T/J = 0.1 (circles), T/J = 0.08 (filled triangles), and 
T/J - 0.0625 (open triangles). Right inset: Data for 5 - 0.9 
with contributions stemming from the lower and upper branch 
shown by dashed-dotted and dashed lines, respectively. 



Examples for the off-diagonal element of the dynamic 
SSF matrix are shown in Fig. [T2] As in the case of 
S(q, q, lu) the two-magnon continuum determines the re- 
gion where S(q, q + n, uj) ^ with edge singularities at 
the boundaries. Moreover we find magnon peaks of dif- 
ferent signs following the nragnon dispersions. 



V. STATIC SPIN-STRUCTURE FACTOR 

The static spin structure factor (static SSF) matrix 
can be easily obtained from the dynamic SSF via 



oo 

S{q 1 q')= I —S(q,q',u). (33) 



We first want to discuss the diagonal elements of this 
matrix. We divide S(q,q) = S~(q,q) + S + (q,q) with 
S~(q,q) (S + (q,q)) being the part which stems from the 
lower (upper) branch of the dispersion. These contribu- 
tions are given by Eqs. (|28[) with a = a' and a ^ a' (see 
also Fig. Ill| . Performing the frequency integrals we find 



S (q,q) 



2N 



o-e{±} k 

x [l + n B (£+*)]. 



(34) 



and 



S + (q,q) = — J2 £[! + *(*, 8 + *)]n B (tf) 



erG{±} k 

x [l + n B (C^) 



(35) 



S ^ <?) - ^ EJ < n+ q +k + n k n q+k <S>{k,q + k) I + S. 

(36) 
Alternatively, the diagonal elements of the static SSF can 
be obtained directly from 






iq{j-i) 



(37) 



However, within the latter approach the distinction of 
the contributions stemming from the lower and the upper 
branch of the dynamic SSF is less transparent. 

In Fig.[l3]thc diagonal structure factor S(q, q) is shown 
for various dimerizations at T/J = 0.1 with 5 = 1. The 
agreement between the MSWT results and the numerical 
TMRG data is good even for high dimerizations. An in- 
creasing dimerization leads to an increase of the linewidth 
of the central peak at q = 0. At high dimerizations a 
shoulder in S(q,q) appears that extends almost up to 
q = 7T. At low tenrpcratures a second maximum develops 
near q = n (see left inset of Fig. [13]). For decreasing tem- 
peratures the peak shifts to higher momenta leading to 
a sharp maximum close to the boundary of the Brillouin 
zone for very low temperatures and high dimerizations. 

For a dimerized FM chain with 5 < 1 such a high mo- 
mentum structure is expected due to the commensurate 
modulation of the exchange constant.— The naive expec- 
tation of a peak right at it indicating the dimerization of 
correlation functions is, however, not fulfilled. Instead, 
we find a small peak with momentum shifted away from 
q = it. This can be readily understood by considering 
S~(q,q) and S + (q,q) separately as is done in the right 
inset of Fig. [13] S~(q, q) is peaked at q = and S + (q, q) 
reveals a maximum at q — tt. However, at sufficiently 
high dimerizations the S~(q,q) contribution also yields 
a significant contribution for q > n/2. Therefore the 
two contributions add up to a peak which is shifted with 
respect to q = tt. 

One may also wonder why the dimerization does not 
show up more pronounced in the diagonal elements of the 
static spin structure factor matrix, i.e., why the peak at 
q = has much more weight than the high momentum 
peak even at strong dimerizations. To investigate this 
property, we insert Eq. (|14p into Eq. (|37p dividing the 
SCF into a uniform and an alternating part. Performing 
the sum over r to lowest order we end up with 



S(q,q) 



1 f A 



D 



with 



and 



6U 2 + er (i-k) 2 + C s 



U ( 1 



(38) 
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FIG. 14: (Color online) Off-diagonal elements of the static 
spin structure factor matrix for S — 1, T/J = 0.1, and various 
dimerizations as a function of q. The symbols show numerical 
data obtained by TMRG, whereas the solid lines are results 
from MSWT, Eq. flU). 



Thus, as expected, we find two Lorentzians at q = and 
q = 7r, respectively. However, the ratio 



S(ir,ir) v$t s S 2 



5(0,0) l + vgtgd 2 



(39) 



is small for effective temperatures fjCl. For large tg, on 
the other hand, we have £g -C 1 and the two Lorentzians 
become very broad. 

Finally, from the left inset of Fig. [13] we see that the 
peak position of the maximum at high momentum is tem- 
perature dependent. On the other hand, we find that for 
ts and 5 fixed the peak position does not depend on S. 
To understand this puzzling feature we derive an ana- 
lytical expression for the position of the high momen- 
tum peak. Keeping only contributions from «b(CaT) m 
Eq. (|31)|) . expanding $(g, k) around q — ir and k = 0, 
and performing a saddle-point integration we can eval- 
uate the remaining expression. Using Eq. (|16p we find 
that the high-momentum peak is located at 



(40) 



which explains why the peak position does not depend 
explicitly on S. 

Next, we turn to the off-diagonal elements of the static 
SSF matrix. From integrating Eq. (f3l"|) over frequency we 
find 




S(q,q + 7r) = — = $^tf(fc,g + *> 



-/.■■ 



(41) 



The off-diagonal elements vanish in the undimerized case 
and are finite in the dimerized case. Thus they provide 
a rigorous criterion for dimerized spin correlation func- 
tions. In Fig. Q3] this expression is shown for 5 = 1, 



T/J = 0.1, and various dimerizations. S(q,q + n) is an 
odd function of momentum transfer q, and is a symmetric 
function with respect to q — n/2. It exhibits a two peak 
structure which is more pronounced the higher the dimer- 
ization is. In complete analogy to the high-momentum 
peak of 5(<7, q) we can again determine the positions of 
the maxima of S(q, q + 7r) and find 



± 



1 
Is 



(42) 



Hence, as in the previous case, the peak positions of 
S(q, q + tt) are determined by the spin quantum num- 
ber 5 and the reduced temperature ts only. In particluar 
we observe that the distance from the center (boundary) 
of the Brillouin zone of the low (high) momentum peak 
is given by the inverse correlation length £7 . 

To summarize, dimerization is signaled by the appear- 
ance of finite intensity in the off-diagonal spin structure 
factor. On the other hand, the conventional (diagonal) 
structure factor reflects the FM spin-Peierls dimeriza- 
tion mainly by quantitative changes, in particular by 
an increase of the lincwidth of the q = peak. As a 
qualitative change, a high-momentum maximum at large 
dimerizations and low-temperatures develops in the diag- 
onal structure factor, however, this structure always has 
a very small intensity. 



VI. CONCLUSIONS 

We have shown that an alternation of the magnetic ex- 
change can also occur in ferromagnetic spin chains. Pos- 
sible mechanisms are a coupling to orbital fluctuations, 
as in YVO ^ 22 ' 23 ^ 8 , or to lattice degrees of freedom. In 
both cases the dimerization has to be thermally activated 
leading to a gain in magnetic energy ~ — T 3 ' 2 <5 2 . In an 
adiabatic approximation the loss in potential energy is 
~ KS 2 where K is an effective elastic constant. A dimer- 
ized phase due to coupling to the lattice is therefore only 
possible as a finite-temperature phase and only if K is 
smaller than some threshold value. 

As an analytical method to calculate the thermody- 
namic properties of a dimerized ferromagnetic Heisenberg 
chain we have used a modified spin- wave theory. By com- 
paring with numerical data obtained by the transfer ma- 
trix renormalization group, we have shown — as a central 
result of this paper — that the modified spin- wave theory 
yields excellent results for a wide range of dimerization 
parameters and temperatures. 

We found that both the spin and the dimcr correla- 
tion functions decay exponentially at finite temperatures 
with the same correlation length £g = c(T, 6)JS 2 /T. The 
correlation length £s decreases with increasing dimeriza- 
tion 6. Next-leading corrections to the 1/T behavior have 
been calculated explicitly for spin 5 = 1/2 and 5=1 and 
even these corrections (of order 1) described by the coef- 
ficient c(T, S) have been shown to be in excellent agree- 
ment with the numerical data. For dimerizations 5 < 1 
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the magnetic susceptibility diverges as \ ~ (1 — $ 2 )/T 2 , 
i.e., the susceptibility at low temperatures is suppressed 
with increasing dimerization. For S — 1 we have a system 
of decoupled dimers and x ~ 1/T". 

Similarly to X-ray scattering experiments, where the 
incident photon is given by a coherent superposition of 
the forward diffracted and the Bragg reflected wave, we 
have formulated the dynamic spin stucture factor for 
neutron scattering experiments as a matrix. The diag- 
onal elements correspond to the standard dynamic spin 
structure factor while the off-diagonal elements include 
a reciprocal lattice vector. The dimerization leads to 
two magnetic excitation branches which are clearly visi- 
ble in the diagonal elements of the dynamic spin struc- 
ture factor matrix. Within the spin-wave approximation 
used here, this quantity is zero outside the two-magnon 
continuum. At the boundary of the two-magnon con- 
tinuum edge singularities are formed due to a diverg- 
ing density of states. These singularities might become 
smooth once higher magnon excitations are included but 
we expect that a peak-structure remains along the two- 
magnon boundary. Experimentally, the edge singulari- 
ties can only be resolved if they carry sufficient spectral 
weight. This is only the case if they occur close to a 
magnon peak. For the off-diagonal elements of the dy- 
namic spin structure factor matrix we also find magnon 
peaks with different signs for the upper and lower magnon 
branch. The two-magnon continuum determines the re- 
gion where the off-diagonal elements of the dynamic spin 
structure factor is nonzero with edge singularities occur- 
ing at the boundaries. 

The diagonal elements of the static spin structure fac- 



tor matrix exhibit q = peaks with Lorcntzian line- 
shape. The linewidth of the latter increases with increas- 
ing dimerization. Furthermore, for strong dimerizations 
a high momentum peak is observed. We find that the po- 
sition of the peak at high momenta depends solely on S 
and the reduced temperature t$ = T/(JS(1 — S 2 )). The 
dimerization leads to a doubling of the unit cell which 
manifests itself by the appearance of a finite off-diagonal 
structure factor. The off-diagonal elements of the static 
spin structure factor matrix show a two-peak structure 
which becomes more pronounced the higher the dimer- 
ization is. The distance of the maxima with respect to 
the center and boundary of the Brillouin zone are given 
by the inverse correlation length. 

To conclude, a spin-Peierls state can be realized in fer- 
romagnetic spin chains but only at finite temperatures. 
Most promising candidates for this new state are, as we 
have argued, monatomic chains on surfaces where the ef- 
fective elastic constants arc expected to be rather small. 
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